Onset of superconductivity in a voltage-biased NSN microbridge 
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We study the stability of the normal state in a mesoscopic NSN junction biased by a constant 
voltage V with respect to the formation of the superconducting order. Using the linearized time- 
dependent Ginzburg-Landau equation, we obtain the temperature dependence of the instability line, 
Vinst(T), where nucleation of superconductivity takes place. For sufficiently low biases, a stationary 
symmetric superconducting state emerges below the instability line. For higher biases, the normal 
phase is destroyed by the formation of a non-stationary bimodal state with two superconducting 
nuclei localized near the opposite terminals. The low-temperature and large- voltage behavior of the 
instability line is highly sensitive to the details of the inelastic relaxation mechanism in the wire. 
Therefore, experimental studies of VSnst(T) in NSN junctions may be used as an effective tool to 
access parameters of the inelastic relaxation in the normal state. 

PACS numbers: 74.40. Gh, 74.78. Na, 72.15.Lh, 72.10.Di 



Nonequilibrium superconductivity has being attract- 
ing significant experimental and theoretical attention 
over decades [IHll, ranging from vortex dynamics Q to 
the physics of the resistive state in current-carrying su- 
perconductors ■ It was recognized long ago that 
a superconducting wire typically has a hysteretic current 
voltage characteristic specified by several "critical" cur- 
rents. In an up-sweep, a current exceeding the thermody- 
namic dcpairing current, I C (T), docs not completely de- 
stroy superconductivity but drives the wire into a nonsta- 
tionary resistive state [111 ] , with the excess phase winding 
relaxing through the formation of phase slips [12]. The 
resistive state continues until I2{T) > I C (T), when the 
wire eventually becomes normal. In the down-sweep of 
the current voltage characteristic, the wire remains nor- 
mal until I\(T) < l2(T) when an emerging order param- 
eter leads to the reduction of the wire resistance. 

The theoretical description of a nonequilibrium super- 
conducting state is a sophisticated problem, requiring 
a simultaneous account of the nonlinear order parame- 
ter dynamics and quasiparticle relaxation under nonsta- 
tionary conditions. The resulting set of equations is ex- 
treme ly comp licated 0, 0] and can be treated only numer- 
ically 13l4l5j (even then the stationarity of the supercon- 
ducting state is often assumed HI, 14 1). A more intuitive 
but somewhat oversimplified approach is based on the a 
time-dependent Ginzburg-Landau (TDGL) equation for 
the order parameter field A(r, t). The TDGL approach 
which is generally inapplicable in the gapped phase [l6| , 
can be justified only in a very narrow vicinity of the crit- 
ical temperature, T c , provided that the electron-phonon 
(e-ph) interaction is sufficiently strong to thermalize 
quasiparticles [l7|. These generalized TDGL equations 
are analyzed numerically in Refs. (Bl [l8j. 

While the applicability of the TDGL equation in the 
superconducting region is a controversal issue, its lin- 



earized form can be safely employed to find the line 
Iinst(T) of the absolute instability of the normal state 
with respect to the appearance of an infinitesimally small 
order parameter A(r, i) [l(| EH, If the transition to 
the superconducting state is second order, then I\(T) 
coincides with L lns t(T). Otherwise the actual instability 
takes place at a larger Ii(T) > Ii ns t{T). In both cases, 

-^inst 

(T) gives the lower bound for h(T). 

Previous results 1(], HI for the instability line of a su- 
perconducting wire connected to normal reservoirs (NSN 
microbridge) have been obtained in the limit of quasi- 
cquilibrium, when strong e-ph relaxation renders the dis- 
tribution function locally thermal. This approximation 
breaks down for low-T c superconducting wires shorter 
than the e-ph relaxation length, l e -ph(T c ) (e.g., for alu- 
minium, l e - P h(T c ) w 40 (im [21]). Such systems have re- 
cently been experimentally studied in Refs. [IH (Al) and 
p2l [23| (Zn; reservoirs may be driven normal by a mag- 
netic field). It was found that for sufficiently large biases 
superconductivity arises near the terminals throu gh a 
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second-order phase transition, with Ii(T) = Ijnst(T') 

In this Letter we study the normal state instability line 
in an NSN microbridge biased by a DC voltage V, relax- 
ing the assumption of strong thcrmalization. For small 
biases, eV <C T c , the instabilit y li ne is universal and we 
reproduce the results of Refs. jloL The universality 
breaks down for larger biases, where we obtain Vj ns t(T) 
as a functional of the normal state distribution function 
and analyze it for various types of inelastic interactions. 

We model the NSN microbridge as a diffusive wire of 
length L coupled at x — ±L/2 to large normal reservoirs 
via transparent interfaces. The terminals are biased by 
a constant voltage V. The wire length, L, is assumed 
to be larger than the zero-temperature coherence length, 
(0 = \/nD/8T c o, where D is diffusion coefficient, and T c o 
is the critical temperature of the infinite wire. The equi- 
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librium critical temperature, T c — T c o(l — 7r 2 £ 2 /L 2 ), is 
smaller than T c q due to the finite size effect [24[ . In what 
follows we neglect superconducting fluctuations pjf . 

General stability criterion. — An arbitrary nonequilib- 
rium normal state becomes absolutely unstable with re- 
spect to superconducting fluctuations if an infinitesimally 
small order parameter, A(r, t), does not decay with time 
but evolves to finite values. Evolution of A(r, t) is de- 
scribed by the TDGL equation, and it suffices to keep 
only the linear term to judge its stability. In dirty super- 
conductors, the linearized TDGL equation can be read- 
ily derived from the Keldysh u-model formalism (26l - [28| 
or dynamic Usadel equations [3] by expanding in A. It 
takes the form (Lp)^ 1 * A = 0, where is the in- 

verse fluctuation propagator, and convolution in time and 
space indices is implied. In the frequency representation, 
is an integral operator in real space specified by 
the kernel 
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S - J j L + ^ I dEF(E,r)C u -2 E (r,r'), (1) 



where A is the dimcnsionlcss BCS interaction constant, 
u>d is the Dcbye frequency, and C stands for the retarded 
Cooperon, C e (r,r') = (r|(— DV 2 — ie)~ 1 \r'), vanishing at 
the boundary with the terminals. 

The operator (jTJ) is a functional of the normal-state 
nonequilibrium electron distribution function, F(E,r), 
which contains information about the mechanism for in- 
elastic relaxation in the wire. The distribution function 
should be determined from the kinetic equation 



DV 2 F(E,r) +l e - e [F] +l e - ph [F] = 0, 



(2) 



with X e ~ e [F] and I e_ph [F] being the electron-electron 
(e-e) and e-ph collision integrals, respectively. The cor- 
responding energy relaxation lengths, l e . e (T) cx T -1 / 4 
and l e - P h(T) cx T~ 3 / 2 , behave as a negative power of the 



temperature T in quasi-equilibrium [21 



In the absence of inelastic collisions, the kinetic equa- 
tion ([5]) is solved by the "two-step" function [29l . |3(| : 

F(E,x) = (1/2 -x/L)F L {E) + (1/2 + x/L)F R {E). (3) 

The distribution functions in the terminals, Fl^r(E) = 
Fq(E ± eV/2), are given by the equilibrium distribution 
function, F (E) = tanh(£/2T), shifted by ±eV/2 (e > 
0). In the opposite case of strong inelastic relaxation, the 
distribution function takes the form 



F in (E,x) = tanh[(£ - e</>(x))/2T(x)], 



(4) 



where <fi(x) = Vx/L is the potential in the normal state, 
and T(x) is the effective temperature. For strong lattice 
thermalization (l e - p h <C L <C l e -e), T(x) = T. For the 
dominating e-e scattering (Z e _ e <C L -C / e -ph), T 2 (x) = 
T 2 + (3/4tt 2 )[1 - (2x/L) 2 ](eV) 2 0. 

The order parameter evolution governed by the lin- 
earized TDGL operator (JTJ) can be naturally described 
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FIG. 1. (Color online) (a) Real (solid blue line) and imaginary 
(dashed red line) parts of the lowest eigenvalues Ai,2(i>) of the 
Hamiltonian ©. The spectrum is entirely real until v — v c « 
49.25. (b)-(e) Spatial dependence of the absolute values of the 
eigenfunctions, |^ji(x)| (blue) and |^2(x)| (red), for v — 0.8n c , 
1.2 v c , 5v c , and 10 v c respectively. 



in terms of the eigenmodes Afe(r)e _Mfet annihilated by 
(L^) _1 . The normal state is stable provided Imcuk < 
for all eigenmodes. Below we analyze the spectrum of 
the operator (JTJ) and determine the instability line in the 
coordinates T— V. In the general case, the spectrum can 
be obtained only numerically. Analytical considerations 
become possible when the operator ([1]) may be linearized 
in to: (L^)^ 1 = iruj — H. The instability occurs when 
the real part of the lowest eigenvalue of the operator % 
turns to zero. 

Weak-nonequilibrium regime. — In the limit of low bi- 
ases, eV <C T c , the deviation from equilibrium is 
small everywhere in the wire and the distribution func- 
tion acquires a universal form, F(E,x) ~ Fq(E) — 
Fq(E) ecj)(x) , regardless of the relaxation mechanism. 
Then Eq. Q takes the form (L^)" 1 = muj/%T - 
ln{T/T c0 )- (to/ L) 2 H V , with 



-ds + 2ivx, 



x £ 



-1/2,1/2]. 



(5) 



The Hamiltonian % v describes quantum-mechanical mo- 
tion in an imaginary electric field, v = eV/E^b., on the 
interval x = x/L 6 [-1/2,1/2], where E Th = D/L 2 
is the Thouless energy. Hard-wall boundary conditions, 
^(±1/2) = 0, imposed on wave functions correspond to 
the complete suppression of superconductivity at the con- 
tacts with highly-conducting terminals due to the inverse 
proximity effect. 

The Hamiltonian (O has been recently analyzed in 
Ref. It belongs to a class of non-Hermitian Hamil- 
tonians invariant under the combined action of the time- 
reversal, T: f(x) H> f*(x), and parity, V: ,f(x) H> f(—x), 
transformations. The "PT-symmetry of the Hamilto- 
nian (O ensures that its eigenvalues X n (v) are either 
real or form complex-conjugated pairs, providin g a com- 
plex extension of the notion of Hermiticity [3ll |32|. At 
V = 0, the spectrum is non-degenerate: A„(0) = ir 2 n 2 
(n = 1,2,...). It evolves continuously with v and a 
nonzero Im X(v) arises only when the two lowest eigenval- 
ues, Xi(v) and \2(v), coalesce [see Fig.[T](a)]. This hap- 
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pens at v — v c sa 49.25 [19J, indicating the transition to a 
complex- valued spectrum. For v < v c . the ground state 
of ([5]) is PT-symmetric, and hence 1^1(5)1 = IVh( — 
For v > v c , the "PT-symmetry is spontaneously broken 
and there is a pair of states with the lowest ReA(i>): 
= ipi(x) and iPr(x) = ip2(%) = 4'K^x), shifted 
to the left (right) from the midpoint [see Fig. [T] (b)-(e)]. 

Spontaneous breaking of the "PT-symmetry associated 
with the spectral bifurcation at v = v c explains the ap- 
pearance of asymmetric superconducting states observed 
in numerical simulations [33j and recent experiments 14 j . 
The normal-state instability line, V; ns t(T), is specified im- 
plicitly by the relation 

1 - T/T c0 = (£ /£) 2 ReA 1 (eF ins t(T)/.E Th ), (6) 

and exhibits a singular behavior at the critical bias V* 
given by eV* = v c Etii ~ 50£rh (see inset in Fig.[2|). The 
bifurcation of the instability line occurs at the tempera- 
ture T* w T c0 (l - 28.44 £ 2 /L 2 ). For long wires (L > £ ), 
T* is very close to T c . 

The time dependence of the emergent superconduct- 
ing state is determined by ImAi(v). Below the bifur- 
cation threshold, for VS ns t(T) < V*, the system under- 
goes at V = VinstCT") the transition to a stationary su- 
perconducting state, with the superconducting chemi- 
cal potential being the half-sum of the chemical po- 
tentials in the terminals. This state is supercurrcnt- 
carrying, and can withstand a maximum phase wind- 
ing of 7r achieved at the critical bias V*. For larger 
voltages, V msi (T) > V*, two modes, iPl(x) and ipn(x), 
nucleate simultaneously at Vi nst (T). The resulting bi- 
modal superconducting state is non- stationary since the 
left and right modes feel different electrochemical po- 
tentials and their phases rotate with opposite frequen- 
cies, fii.nOO = tEtil Im Ai (eV/T) . This will result in 
the Josephson generation with the differential frequency 
Q(V) = n R (V) - n L (V) oc y/V-V, as V inst (T) -> V*. 
Though these supercurrent oscillations are locked to the 
superconducting part they may excite oscillations of the 
normal current in the whole circuit. Thus the dc biased 
NSN microbridge may act as a voltage-tunable generator 
of an ac current, with the maximal amplitude of oscilla- 
tions expected in the coherent regime, V — ~ V*. The 
possibility of experimental observation of such a genera- 
tion remains an open problem. 

Incoherent regime. — As the voltage is increased far 
above the bifurcation threshold, Vj ns t(T) 3> V*, the 
cigenmodes ipL,R.{x) gradually localize near the corre- 
sponding terminals, with their size, a(V), becoming much 
smaller then the wire length [see Fig. [T] (b)-(e)]. This is 
the incoherent regime, where the overlap between tpL (x) 
and ipn{x) is exponentially small, supercurrent oscilla- 
tions are suppressed, and nucleation of superconductivity 
near each terminal can be described independently fl4| . 

Using a(V)/L as a small parameter and still work- 
ing in the vicinity of T c , we linearize F (E, x) near the 




FIG. 2. (Color online) Instability voltage as a function of tem- 
perature, Vinst(r), obtained numerically for a wire of length 
L = 15£o for three limiting types of the distribution func- 
tion: without inelastic relaxation (solid blue line), and with 
dominant e-e (dot-dashed line) or e-ph (dashed line) relax- 
ation. Dotted curve illustrates the suppression of V^^ e ^(T) 
by a finite terminal resistance, fi = 0.1 (see text). Inset 
shows the behavior in the vicinity of the bifurcation point 
eV, « 50-Erh « 0.57 T c0 . 



left terminal and reduce Eq. ([T]) to the form: (L^) 1 = 
iir(oj + eV)/8T — \n(T/T c o) — W a , where the operator 

H a = -tld 2 XL + ax L , x L > 0, (7) 

acts on the semiaxis xl = x + L/2 > with the bound- 
ary condition ip(0) = 0. The complex parameter a is a 
functional of the distribution function: 



dE d x F(E - eV/2, x)\ x =_ L/2+a{v) 



(8) 



2{E-iO) 

The ground state of the Hamiltonian Q has the energy 
7o(a£o) 2 ^ 3 and the wave function iPo(xl) = Ai(xi/a — 
7o)/\/a Ai'(— 70), where —70 ~ —2.34 is the first zero 
of the Airy function, and a — (Co/^) 1 ^ 3 i s the nucleus 
size [lo| . For the instability line we get: 

1 - T/T ca = 7o ^ /3 Rea 2 / 3 (eV in8t (T)/T). (9) 

The left and right unstable states rotate with the fre- 
quencies fl L>R (V) = T[eV - CLiiV)}, where Q^V) = 
(8T' c o/7r)7o^ 2//3 Ima 2 / 3 (eU/T) is a small correction to the 
Josephson frequency determined by the electrochemical 
potential of the corresponding terminal. At the instabil- 
ity line, V = Vi aB t{T), the size a of the unstable mode 
is of the order of the temperature-dependent supercon- 
ducting coherence length f (T) - (1 - T/T c0 )- 1/2 £o- 

For long wires (L £0)1 the incoherent regime partly 
overlaps with the weak-noncquilibrium regime. Then for 
eV* <C el^ ns t(T) <C T c , Eq. ((9]) gives a universal answer 



eV inst (T) _ 2 7 / 2 L ( T c q — T 



T, 



<:■() 



£0 V To^cO 



3/2 



(10) 
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which could have also been deduced from Eq. ([6]) us- 
ing the quasiclassical approximation for Xi(v) at d > 1. 
Equation ([§]) exactly coincides with the result of Ref. , 
predicting that superconductivity nucleates near the ter- 
minals at a finite current Ii n st(T) w 0.356 I C {T). 

The position of the instability line in the incoherent 
regime at large biases, eVi nst (T) ^> T c , depends on the 
relation between the inelastic lengths l e . e and l e -ph, the 
wire length L, and the nucleus size a(V). The presence 
of the latter scale, which probes the distribution function 
near the boundaries of the wire, leads to a rich variety of 
regimes realized at different temperatures. 

For the three limiting distributions [Eqs. ((3]) and 
(0J)], the function a(u) can be found analytically: (i) 
a free (it) = [-0(1/2 + iu/2n) - i/>(l/2)]/L for the non- 
interacting case, L <C l e - e ,le-ph, where ip( x ) is the 
digamma function; (ii) a e _ p h(w) = iiru/AL for strong lat- 
tice thermalization, Z e _ p h -C ct(V) «L« / e _ e ; and (hi) 
a e . e (u) = [iitu/A + 3u 2 /2ir 2 ]/ L for the dominant e-e in- 
teraction, Z e _ e <C a(V) <C L <C i e -ph- In case (ii), the 

instability line V^ t h (T) is S ivcn h V Ec l- CHI)- In tlic 
vicinity of T c , the instability lines in cases (i) and (hi) 
are given by: 



eV 



(free) 



(T) 



l cO 




(11) 

(12) 



Countcrintuitively, in cases (i) and (iii) the instability 
current L ms t(T) oc V^ ns t(T)/L has a nontrivial depen- 
dence on the system size, as opposed to Eq. (fT0|) . Such 
a behavior is a consequence of strong nonequilibrium in 
the wire. The limiting curves v£™\t), V^ 11 ' (T), and 
Vj^j. (T) for all temperatures obtained numerically from 
Eq. ((XJ) for the wire with L/£q = 15 are shown in Fig. [2] 
The universal behavior at small biases can be easily seen 
(inset). Since the ratio £/£o is not very large, the insta- 
bility line becomes strongly dependent on the distribu- 
tion function already for V > V* . 

The most exciting feature of our results is the expo- 
nential growth of Vi ns t(T) with decreasing temperature 
in the non- interacting case, Eq. pip . Hence, even a small 
deviation of the distribution function from the two-step 
form J3]) will drastically modify Vj ns t(T). As an example, 
consider the effect of a finite resistance of the normal ter- 
minals. Then the function Fl(E) in Eq. will be re- 
placed by F L (E) = ftF (E+eV/2) + (l-ft)F (E-eV/2), 
where V is the voltage applied to the NSN microbridge, 
and ft = Rt/(Rn + 2Rt) [Rt and Rn are the resistances 
of the N and S part of the junction, respectively]. The 
resulting Vj ns t(T) for ft = 0.1 is shown by the dotted blue 
line in Fig. [5] While Vi ns t(T) is unchanged for small bi- 
ases, it is strongly suppressed compared to (T) for 
large biases. 



Low-temperature behavior. — The exponential growth 
°1 ^inst (T) in the non- interacting case formally implies 
that superconductivity at T = might persist up to ex- 
ponentially large voltages, ln(eV;nst(0)/T , c o) ~ £/£o 3> 1- 
This conclusion is wrong, since inelastic relaxation and 
heating become important with increasing V, even if they 
were negligible at V — 0. To study the low-T part of the 
instability line, we consider here a model of the e-ph 
interaction (e-e relaxation neglected) when the phonon 
temperature is assumed to coincide with the base tem- 
perature of the terminals and e-ph relaxation is weak at 
Tc- le-ph(T c ) ^> L (as in the experiment [14|). 

With decreasing T below T c , the instability line first 
follows Eq. (ITT1) . At the same time, Z e -ph decreases 
and eventually the distribution function in the mid- 
dle of the wire becomes nearly thermal with the ef- 
fective temperature T g. This happens when T g ob- 
tained from the heat balance equation [2l|, (eV/L) 2 ~ 
^eff/ T c'e-ph( T c). becomes so large that l e - P h(T e $) ~ L. 
The corresponding voltage, Vph, can be estimated as 
eVph/T c ~ [Z e _ p h(T c )/L] 2 / 3 . Consequently, the exponen- 
tial growth (|TT|) persists for voltages V* < V < V p b, 
corresponding to the temperature range T p h < T < T*, 
where with logarithmic accuracy 1— T p h/T c o ~ {£,q/L) 2 / 3 . 

For higher biases, V > Vph, electrons in the central 
part of the wire have the temperature T e fj. However, 
the parameter a, Eq. ©, is determined by the distribu- 
tion function in the vicinity of the terminals which is not 
thermal. Matching solution of the collisionless kinetic 
equation for < xl < le-ph(T e ff) at the effective right 
"boundary", xl = le-ph(T e s), with the function with 
T(x) = Toff, we obtain a ~ l/Z e -ph(T e g). Therefore, for 
V ^ Vph we get with logarithmic accuracy: 



eMnst(T) L 



T, 



^e-ph(-^c) 



2/3 



T c0 ~T 
T, 



5/2 



(13) 



Equation (|T3|) corresponding to the case a(V) <C l e -ph "C 
L is different from the expression (| 10[) when phonons are 
important already at T c , and l e - P h <C a(V) <C L. The 
scaling dependence of Eq. (fT3)) on L indicates that the 
stability of the normal state is controlled by the applied 
current, similar to Eq. (|10p. At zero temperature the 
instability current exceeds the thermodynamic depairing 
current by the factor of {l e -ph(T c )/£o) 2/3 > 1. 

Discussion. — Our general procedure locates the abso- 
lute instability line, Vi ns t(T), of the normal state for a 
voltage-biased NSN microbridge. Following experimen- 
tal data [3] we assumed that the onset of superconduc- 
tivity is of the second order. While non-linear terms in 
the TDGL equation are required to determine the order 
of the phase transition [34[ , we note that were it of the 
first order, its position would be shifted to voltages higher 
than V inst (T). 

In the vicinity of T c , the problem of finding Vi ns t(T) can 
be mapped onto a one-dimensional quantum mechanics 
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in some potential U(x). For small biases, eV <C T c q, 
the potential U(x) does not depend on the distribution 
function details, explaining universality of the instability 
line, including the bifurcation from the single-mode to 
the bimodal superconducting state at eV ~ 50i/rh [H| 
and nucleation of superconductivity in the vicinity of the 
terminals for larger biases [Toj ] . 

For eV > T c q, the potential U(x) becomes a func- 
tional of the normal-state distribution function, produc- 
ing Vi ns t(T) that is strongly sensitive to inelastic relax- 
ation mechanisms in the wire. For the dominant e-ph 
interaction, the instability is controlled by the electric 
field £ = V/L [Eqs. (JTUJ) and ((15)1 ]. while in the opposite 
case [Eqs. (JTTJ) and 1(12)1]. the instability cannot be solely 
interpreted as current or voltage-driven. At zero temper- 
ature, the (nonuniform) superconducting state can with- 
stand a current which is parametrically larger than the 
thermodynamic depairing current. 

High sensitivity of Vi ns t(T) to the details of the dis- 
tribution function opens avenues for its use as a probe 
of inelastic relaxation in the normal state. The shape of 
Vi ns t(T) can be further used to determine the dominat- 
ing relaxation mechanism and extract the corresponding 
inelastic scattering rate. 
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J. C. W. Song, and D. Y. Vodolazov for discussions. 
This work was partially supported by the Russian Federal 
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